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Abstract 

In this paper we present a non-singular black hole model as a possible 
end-product of gravitational collapse. The depicted spacetime which is 
type [II, (II)], by Petrov classification, is an exact solution of the Ein- 
stein equations with and contains two horizons. The equation of state 
p r (p) in the radial direction, is a well-behaved function of the density 
p (r) and smoothly reproduces vacuum-like behavior near r = while 
tending to a polytrope at larger r, low p values. The final equilibrium 
configuration comprises a de Sitter-like inner core surrounded by a family 
of 2-surfaces E of matter fields with variable equation of state. The fields 
are all concentrated in the vicinity of the radial center r — 0. The solution 
depicts a spacetime that is asymptotically Schwarzschild at large r, while 
it becomes de Sitter-like as r — > 0. Possible physical interpretations of the 
macro-state of the black hole interior in the model are offered. We find 
that the possible state admits two equally viable interpretations, namely 
either a quintessential intermediary region or a phase transition in which 
a two-fluid system is in both dynamic and thermodynamic equilibrium. 
We estimate the ratio of pure matter present to the total energy and in 
both cases find it to be virtually the same, being ~ 0.83. Finally, the well- 
behaved dependence of the density and pressure on the radial coordinate 
provides some insight on dealing with the information loss paradox. 

Keywords: non singular black holes 

* NAS/NRC Senior Research Associate; e-mail: mrmsps@rit.edu 



Pacs: 



04.70.Bw, 04.20.Jb 



1 Introduction 

It is generally accepted that the formation of a black hole proceeds with the 
formation of a marginally trapped surface (the apparent horizon) into which 
matter is collapsing and one that encloses a region of space interior to which 
the light cones arc tipped [1] so much so that even outgoing null geodesies 
proceed toward ever decreasing values of the radial coordinate. The famous sin- 
gularity theorems [2] [3] [4] , then posit that under certain assumptions on energy 
conditions of matter, the eventual state of the collapse is an infinite density 
singularity surrounded by an event horizon (defined by the outermost trapped 
surface). It is, however, a fact that the existence of such a singularity is nei- 
ther easy to verify nor refute observationally because of the limits set by the 
black hole event horizon. In the absence of such observational evidence, the 
compelling logic in the arguments leading to the singularity theorems has gen- 
erally shaped the widespread conviction that indeed a black hole must contain 
a physical singularity. 

As is known, however, the picture of a singular black hole opens up sev- 
eral puzzles. For example, the existence of spacetime singularities results in 
irretrievably total information loss [5]. Moreover, the infinite tidal forces that 
result from such a singularity formation lead to a breakdown in the descriptive 
power of General Relativity [6]. As a result it is generally believed that the 
perceived final collapse state implied by the singularity theorems could be but 
a manifestation of the incompleteness of classical gravity and that a correct 
(quantum) theory of gravity will dispense with this singular final state and ad- 
mit a state of finite (albeit very high) curvature. An altogether different view is 
that the singularities implied by the aforementioned theorems simply reflect our 
lack of knowledge of the properties of matter under extreme conditions. Thus 
it has been suggested that the singularity theorems could be circumvented if 
(some of) the conditions imposed on the properties of collapsing matter were 
relaxed. Indeed modern frameworks (like string theory) attempting to formulate 
a quantum theory of gravity seek the existence of a fundamental length scale 
and hence a singularity free spacetime. There are also philosophical questions 
as to whether spacetime can indeed accommodate singularities of any kind, be 
they from gravitational collapse or vacuum- induced [7] . 

As a result of these and other considerations a renewed interest has grown 
lately with regard to spacetime singularities and their existence [8] . Indeed the 
concept of non-singular collapse dates back to Sakharov's consideration of the 
equation of state p = —p for a supcrdense fluid and Gliner's view [9] that such 
a fluid could be the final state of gravitational collapse. Such ideas are now 
attracting a growing amount of attention. For example, Ellis [10] has argued 
that in the case of a closed universe a trapped surface does not necessarily lead 
to a singularity. Recently the present authors have investigated an analogous 
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argument with regard to black holes in [11]. It was pointed out in [11] that under 
physically acceptable energy conditions the expansion of the outgoing mode 
which assumes increasingly negative values in a trapped region can eventually 
turn around to (even) assume a positive signature inside a black hole. Such 
behavior suggests the existence some interior region, surrounding the core, which 
is not trapped. In turn, this feature implies the absence of a singularity in such 
a spacetime. In recent years, several non-singular black hole solutions have been 
found. Dymnikova, for example, [12] has constructed an exact solution of the 
Einstein field equations for a non-singular black hole containing at the core a 
fluid dc-Sittcr-likc fluid with anisotropic pressure. This solution is Petrov Type 
[(II), (II)] by classification 1 . Another line of investigation in this area, initiated 
by Markov [13], suggests a limiting curvature approach. This idea has been 
explored further by several authors, see for example [14] [15] [16] [17] [18] [19]. 

A common feature in all these treatments is that the geometry of the space- 
time in question is Schwarschild at large r as expected by Birkoff's theorem, 
and de Sitter-like at small r values as implied by the p = —p equation of state 
assumed to prevail near r = 0. These approaches can be put in two broad 
classes, (a) Those for which the transition to this "exotic" state of matter can 
be placed well inside the Schwarzschild horizon [12] and (b) those which replace 
the entire volume of the Schwarschild metric interior to r = 2M by a substance 
with the p = —p equation of state [6], dispensing in this way with the presence 
of the horizon altogether and the "information paradox" problems it engenders. 

The problem of transition from the usual matter equation of state to that 
appropriate for these non-singular solutions is usually not addressed at all. This 
problem is particularly acute for models of the second class discussed above, 
since the density at horizon formation scales like p ~ 10 16 (M/Mq) -2 g cm~ 3 
and for objects of mass M > 10 8 M Q ), is not greater than that of water, which is 
known to have an equation of state very different from p = —p. Moreover, issues 
to do with direct matching of an external Schwarschild vacuum to an interior de 
Sitter have previously been discussed [20] based on junction conditions [21]. In 
fact, in terms of gravitational collapse the associated difficulties have been used 
to suggest modifications in some treatments [22] . The present work is motivated 
by these questions and is undertaken as a first effort to provide some answers 
by constructing models using an explicit equation of state with the desired 
properties. It would appear that the junction constraints [21] when applied to 
the end-product of non-singular gravitational collapse imply the matter fields 
across the Schwarschild/de Sitter boundary should have a radially dependent 
equation of state of the form 1 < w(r) < — 1 that smoothly changes the matter- 
energy from a stiff fluid to a cosmological constant. 

The model sketched in this paper describes the geometry of a body that 
passed during its collapse through a stiff fluid state to settle into the final state 
of a non-singular black hole. The solution allows for a matter fluid region (with 
a family of equations of state < w (r) < 1 which envelopes an intermediate 

1 We shall present a solution with the same asymptotic behavior as in [12] which however 
is different in both classification (it is Petrov Type [II, (II)] and physical interpretation. 
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region with a family of equations of state — 1 < w (r) < which in turn envelopes 
a de Sitter like region with w = — 1 at the center. Our model differs both 
from the traditional singular black hole solutions which introduce matter at 
the singularity and from the relatively new non-singular solutions which usually 
introduced only a de sitter-like spacetime inside the black hole. The main feature 
of our model is that it 

(1) introduces gravitating matter inside a non-singular black hole and 

(2) offers a reasonable explanation of how part of this matter can evolve 
towards a de Sitter-like vacuum to provide the radial tension or negative pressure 
that supports the remaining matter fields against the formation of a singularity. 

The rest of the paper is organized as follows. In section 2 we summarize the 
equations to be solved. We also construct the working equation of state and 
discuss some of its desirable features. In section 3 we solve the Einstein Field 
Equations for the desired geometry and obtain an exact solution. In section 4 
we highlight the physical implications of the model. Section 5 concludes the 
discussion. 

2 Problem formulation 

2.1 The field equations 

It is assumed, for simplicity, that the collapsed object can be reasonably de- 
scribed by a spherical, static geometry. The desired line element then takes the 
general form 

ds 2 = -e v ^dt 2 + e x(r Ur 2 + r 2 (dO 2 + sin 2 6d<p 2 ) , (2.1) 
where A (r) and v (r) are to be determined from the Einstein Field Equations 

G^u = SttGT^u. (2-2) 

Under the assumed spacetime symmetry, Eqs.2 reduce to 

^(7-^)+ ^ = 8*<?P(r), (2.3) 

e ~ A (7 + ^) " ^ = ^ G Pr W , (2-4) 

and 

x (v" AV v 12 v' -V\ a 

Here p (r) = Tq is the energy density, p r = T/ is the radial pressure and 
(pg = Tf ) = (p<p = X3 ) = p± is the tangential pressure. In this treatment the 
fluid is anisotropic with T 2 2 = T 3 3 ^ T/ and T/ ^ T ° ^ T%. The spacetime is 
thus Petrov Type [II, (II)], where ( ) implies a degeneracy in the eigenvalues of 
the Weyl tensor. 
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2.2 Equation of state 

In the remaining part of this section we construct the working equation of state 
p r (p) in the radial direction which reproduces the characteristics of the collapsed 
body as highlighted in the previous section. We justify the choice by pointing 
out its desired features. On the other hand, the tangential pressure equation of 
state p± (p) will be derived later from the Einstein Field Equations (Eq. 2.5) 
using the radial pressure function p r (p) to be constructed. 

The main desired features are that (1) the equation of p r (p) change from 
the one associated with the usual matter, at low densities, to that associated 
with the de Sitter geometry at sufficiently high densities near the center of the 
configuration and; (2) that this change be well-behaved. To this end we consider 
an equation of state that takes the general form 2 where to and 1/n are, yet to be 
determined, positive real numbers and a is a free parameter to be constrained. 
This equation implies presence of a 



Pr (p) 



a-(a + l) 



P 



l/r. 



P, (2-6) 



maximum limiting density p max concentrated in a region of size r = A / — , 

y ^Pmax 

the core region of the configuration. One notes that the main desirable fea- 
tures are readily apparent in the functional form of Eq. 2.6. At low densities, 

{ p p J < the equation of state reduces to that of a polytrope of index 

n (p r oc p 1+1 / n ) while for ( j^— ) > the pressure decreases to eventually 

approach that of the vacuum p r = —p max as p — > p max . We look for the simplest 
form of Eq. 2.6 by judiciously constraining both the indices to and n and the 
parameter a. 

Since Eq. 6 already satisfies the general features desired by the model, its 
specific form, based on an appropriate choice of the indices to and n and the 
parameter a is made through the demand that it must satisfy the following 
basic conditions. 

(1) It must have no pathologies, i.e. 

(i) the sound speed can not be a maximum at p = 

(ii) in order to rule out superluminal behavior, 0^1, the maximum 
sound speed, i.e. the value of |2 a t the point where ^§ = must be given by 
& \ - r 2 - r - 1 

dp \Pstiff s 

(2) It must satisfy minimum acceptable energy conditions, namely the Weak 
Energy Condition, p > 0, p+p r > and the Dominant energy Condition, p > 0, 
Pr e [-p + p\. 

The [to = 1, 1/n = 0] and the [m = 2, 1/n = 0] cases are both pathological 
(they do not satisfy (i) above) and will therefore be discarded. The next simplest 
potential choices are [m = l,n = 1] and [to = 2, n = 1]. Both these forms, ap- 
parently, manifest no pathologies of the previous cases. In both cases the sound 

2 This choice may not be unique. 
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speed vanishes c s — ► at vanishing density p — * 0. Furthermore for /? > the 
sound speed c s is initially an increasing function of p, as expected in reality??. 
One can therefore demand in each case (i.e. for fixed m and n values ) that the 
maximum sound speed coincide with that of light i.e. the cxtremum sound speed 
at = be given by ^ = 1 in order to rule out superluminal behavior. These 
conditions also suffice to determine the parameter a and allow the fulfillment 

,2 j 

of the energy conditions in (2). Solving -jS = and -fi = 1 simultaneously for 

the [m — l,n= 1] case, i.e. fixes a at | — or |\/33 + | (we keep only 

the positive solution). Similarly for the [m = 2, n = 1] one fixes a to a = 2. 
213 5. We will select the quartic function corresponding to m = 2, rz = 1 as the 
representative form for our model equation of state, in part because, unlike the 
cubic form, it unambiguously provides only one possible value of a and makes 
the interpretation simpler. Thus Eq. 2.6 now takes the form 



with the constraint that a = 2.213 5. It is the equation of state we adapt 
throughout the remaining part of the paper. 

A few words on the behavior of our assumed equation of state. While 
no explicit microscopic prescription leading to its form is presently given, we 
consider that its softening past the "stiff" (c 2 s = dp/ dp = 1) state and its 
eventual conversion to an equation of state appropriate to that of the vacuum 
(p = —p), is effected by the coupling of matter to a scalar field akin to that 
of Higgs that gives rise to particle masses. The dominance of the pressure 
by terms which render it negative is considered to imply that eventually the 
energy associated with the self interaction of this field provides the dominant 
contribution to the energy momentum tensor, which is assumed to be that of a 
perfect fluid. 

In this paper no attempt is made to follow the dynamical evolution of the 
collapse. Instead it is assumed that the collapsed body has already reached static 
equilibrium from its collapse. Thus, we only investigate the characteristics of 
the various static 2-surfaces as the radial coordinate decreases from the matter 
surface r = R to the body center r = 0. Since our anticipated solution will allow 
matter fields we can, with no loss of generality, take for initial conditions on the 
collapsed body surface R to be p — 0, p r — 0. We therefore expect to have four 
regions in this spacetime which must be satisfied by the expected solution. 
Region I: Schwarzschild vacuum: R < r < oo, p = p = 
Region II: Regular matter fields: r £ < r < R, p > 0,0 < p < p, 
Region III: Quintessential fields: r < r < r e , < p < p max , — Pmax < P < 0, 
Region IV: A— vacuum: < r < tq, p = —p = — p max . 

The solution to Eqs. 2.1-2.5, must satisfy the following asymptotic condi- 
tions at large r and small r, respectively. 

(i) The spacetime must be asymptotically Schwarzschild for large r, i.e. for 




(2.7) 



G 



R < r < oo 



ds 2 = -[ l- — \dt 2 + - ^Mldr 2 - r 2 (dO + sin 2 0d<p 2 ) . (2.8) 



where for the black hole R is some hypersurface such that R < 2M, M — 
47T J °° p (r) r 2 dr being the total mass. 

(ii) The spacetime must be asymptotically de Sitter 



r 2 



r, 



ds z = - 1 - -j \dt z + t v dr z + r A (dO + sin z 6dip z ) . (2.9) 



('-*) 



for < r < (r < i?). Here, r = y S7rGp — = y X signals the onset of de Sitter 

behavior, where p max = P \r >o is the upper-bound on the density of the fields 

of order of Planck density ppi. 

The asymptotic conditions in Eqs. 2.8 and 2.9 imply there is an interior 
region which includes the family of surfaces E = {S// U S///} and which in- 
terfaces with region I on the outer side and region IV on the inner side. The 
entire spacetime must therefore satisfy regularity conditions at the two inter- 
faces r = R and r = ro. Put simply, such conditions guarantee (i) continuity 
of the mass function and (ii) continuity of the pressure across the interfacing 
hyperfaces. Thus at each interface, i.e. (I, II) and (///, IV) we must have 

[p + -p1, r=ri = o, 

(2.10) 

where r» = {R, r } and +, — refer to the exterior and interior values respectively 
and p — {p r ,Pj_}- A desirable feature of the model (as we find later) is the 
smooth continuity of the density and pressure between region II and region III. 
This feature removes the problem having to match the two regions through 
junction conditions between the regular matter fields and the vacuum-like field, 
since here such conditions will be satisfied naturally. 



3 The solution 

We now solve the Einstein equations 2.3-2.5 for the spacetime of the model. 
Integration of Eq. 2.3 gives 

e- x = 1 - 8n- f p (r') r l2 dr' = 1 - (3.1) 
r Jo r 

where, as stated before, m (r) is the mass enclosed by the a 2-sphere of radius 
r. Further, using Eqs. 2.7 and 11, integration of the Eq. 2.4, gives 

2m(r)\ J M°-(°+P(^) 2 1 (^Vt^ttt 
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where a = 2.213 5. 



Eqs. 3.1 and 3.2 into Eq. 2.1 give 



(l - ^) e r ^dt 2 + r 4^ T ^ 2 + r 2 (d6 2 + sin 2 9d<p 2 ) , (3.3) 



where 

r 



M = / 



8tt 



a - (a+1) 



r — 2m (r) 



pdr (3.4) 



Eq. 3.3 is the general solution for the geometry of our model. It describes 
the spacetime of a non-singular black hole with both matter and a de Sitter 
core. We shall describe below a particular solution resulting from the choice of 
a prescribed density function p (r) . 

To formally complete the solution one must also solve equation 5 for the 
tangential pressure p± = {pe,Pip}- On doing this one finds (see also [23] [24]) 
that 



1 



P±=Pr + nPr + 77 (Pr + P) 



Gm (r) + AirGr^p,, 



r - 2Gm (r) 



(3.5) 



where p r is given by Eq. 2.7. Note that in this model the last term in Eq. 3.6 is 
not vanishing, in general, as is the case considered in some previous treatments 
where it was assumed that p r = —p for all p (see for example [25]). Eq.3.6 is a 
generalization of the Tolman-Oppenhcimcr-Volkoff equation [26]. 

One can now discuss the results of Eqs. 3. 1-3.4. The mass m (r) enclosed a 2- 
surface at any radial coordinate < r < R is given by m (r) — 4ir p (r') r' 2 dr' . 
The entire mass of the body M is given by M — L°° m (r) dr = J °° p (r) r 2 dr. 
It is then clear that outside the mass of the body Tq = 0. Further, as p vanishes 
outside the body, so does the pressure (see Eqs. 2.7 and 3.3). In particular, 
from Eq. 3.3 one observes that in the limit p — > and m (r) — ► M the lapse 
function and the shift vector give e v = e~ A = 1 — 2G r M . Thus outside the mass 
(i.e. for large r) the spacetime becomes asymptotically Schwarzschild, with the 
metric given by Eq. 2.8. Further, in the limit r — > ro, p —> p m ax (and as both 
Eqs. 2.7 and 3.4 show) p r — > -/5 max and p± — > — p max - It follows, therefore 
that for < r < r the fluid has constant positive density p max and constant 
negative pressure — p m ax and takes on the character of a de Sitter spacetime, 
through the equation of state p = -p max - 

In order to establish the actual asymptotic behavior of the spacetime at 
small r one has to choose a density profile for the matter. We choose to 
use the one employed by Dymnikova [12] because of its simplicity and con- 
venient form for integration over the source volume p = p max exp 
The mass m (r) enclosed by a 2-sphere at the radial coordinate 



by m (r) = f Q r ^^{exp 



^0? 



}r' 2 dr' = M 



1 — exp 



giro) 2 

r is then given 
where M 



8 



is the total mass M 



rOO 

Jo P* 



t {exp 



}r 2 dr. The entire mass is essen- 



which in general is close 



o' g 



,1/3. 



tially concentrated in a region of size R ~ (r, 
to r = 0, with its precise value depending on the value of the maximum density 
assumed p max . 

' in Eq. 3.3 



0, it is seen that 



The introduction of the density function p = p max exp 

provides a particular solution for the model. In the limit 
the solution Eq. 3.3 asymptotically leads to Eq. 2.9 as the spacetime becomes 
asymptotically de Sitter. 



4 Physical Interpretation 

Eqs. 3.3, 3.4, 3.5 along with the prescribed density profile, form the solution 
in this model describing the spacetime of a non-singular black hole containing 
both matter and a de Sitter core. A heuristic interpretation of the form of the 
equation of state used has been given in the previous section. The field whose 
presence leads to the asymptotic relation p = —p has not been specified but 
it must be similar to the quintessence field used in cosmology. In this specific 
case, it could in fact be the Higgs field that gives rise to the particle masses. 
As such the maximum density could be roughly p max — yxih-i the latter being 
the Higgs field expectation value. This could be either of order 1 TeV or even 
of order of the Plank mass if, as argued, that should be the typical mass of all 
scalars coupled to gravity. 

The depicted scenario allows two viable, albeit, different interpretations in 
describing the internal structure of such a black hole. The two interpretations 
have the common feature that in both cases a fluid with a negative pressure, in 
the interior, supports the matter fields in the outer region against any further 
collapse. The effect in both cases is to create a non-singular black hole. 



4.1 The quintessential picture 

The first interpretation is directly based on the analytical functional form of the 
equation of state depicted in Eq. 2.7. Here, as one moves from the surface R of 
the collapsed body towards the center, one first encounters region II where from 
the functional form of the equation of state one crosses a family S// of 2-surfaces 
with matter- like w > 0. At the critical point p c in the p — p r space ^ = 
thereafter the pressure p r begins to drop with increasing density p. At some 

point p c = (^ j^pr^ 2 Pmax, the fluid temporarily becomes pressureless w — > 0, 

implying that = G\ — > 0. Beyond this, for smaller values of r one enters 
region III in which one is now crossing a family of 2-surfaces Ejjj for which the 
equation of state is given by —1 < w < 0, smoothly decreasing from w = 
towards w = — 1. In this picture we will refer to this as the quintessential region. 
Evidently exhausts the entire family of quintessential fluids — 1 < w < 0. 
This fluid family provides some of the negative pressure against the implosion 
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of the matter fields in the outer region. The rest of the negative pressure is 
provided by the constant density p max , constant pressure p = -p max , inner core 
< r < r that mimics the cosmological constant. This interpretation is valid 
provided the equation of state p r (p) is, indeed in reality a well-behaved function 
of r for p c < p < p max . 

In this quintessential picture, one can estimate the amount remaining pure 
matter (the total field for which w > 0) as a fraction of the total energy of the 
system. This matter lies in the outer region r e < r < R (i.e. between points B 
and C in Figure 1). Since B is given by, ^ \ Pr — and C is given by p r \ Pe = 0, 
one can infer from Eq. 7 that in this region the density function p is bounded by 



Pc — ( 3 ( Q +i)) Pmax 



<P< 



((«+!)) 



To determine the matter 



mass Mmatter enclosed we use the density profile p (r) = p m ax exp — and 

L r o r g _ 

perform the integral J pr 2 dr with the appropriate limits. It is useful to change 

and we 



the independent variable from r to p. Then r 2 dr = ^r^Tg y^y 
have that Mmatter = 47r / pr 2 dr = f 7rrgr g f£ £ dp. Thus 



) d ( ) 



M 



matter 



a 



a + 1 



(4.1) 



Using a = 2. 213 5 as adopted earlier for our model, gives M mat ter — 0.829 95 (^r 2 r g ) p n 



On the other hand the total mass M — Att f pexp -5 — r 2 dr 
Thus the fractional content of matter in the quintessential picture for the black 



4 2 

2 TTVq Tg /?max • 



hole is ^ 



M 



0.829 95 - 0.83. 



4.2 A two-fluid system 

One can attempt yet a seemingly equally viable interpretation of the black hole's 
internal structure. If one holds as in [6] the view of a phase transition of the fluid 
from the matter-like form to the de Sitter state given by p = — p max , then the 
successive 2-surfaces in the family r < r < r c will represent a two fluid system 
of matter fields and the A— vacuum which in both dynamic and thermodynamic 
equilibrium. The fluid begins as pure matter on the 2-surface r c and gets richer 
and richer in the A— vacuum state as one moves deeper and deeper into the 
black hole. Then the density p r (r) at any such point on any such 2-surface 
shell r < r < R is really a sum of two partial densities 

p = /?(r)p c -[l-/?(r)]/w, (4-2) 

where p c is the critical density just before collapse when p = p c and < [3 < 1. 
The function (3 (r) is the fractional content of pure matter fields to the total 
surface energy contained in an elementary 2-surface shell 4irr 2 dr at ro < r < R. 
For decreasing r, (3 (r) is a decreasing function, evolving from (3 (R) = 1 to 
(3 {ro) — 0. In terms of the partial pressures of the matter fields and the de 
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Sitter vacuum, the pressure at any position r < r < R is then given by 



Pr (r) = 



n=0 



[I3(r)p c -(l-p(r))p n 



(4.2) 



where the o's are positive constants given respectively by 
ao = (1-/3) [a - (a + 1) (1 - (3) 2] 
ai = f3[a + (a + 1) (1-/3)], 
a 2 = /3 2 (a + 1) (1 - /?) , 
a 3 = (3 3 (a + 1) . 



(4.4) 



This expression based on partial pressures will then replace Eq. 2.7 as the 
equation of state p r (p) in the event the evolution of matter fields to a de Sitter 
state involves a spontaneous phase transition. One notes that the pressure p r 
still has the desired asymptotic behavior. Thus in the limit (3 (r) — ► 1, p — > p c = 

a - (a + 1) (-^-) i^-) Pc while in the limit (3 (r) -> 0, p -> -p max as 

expected. 

One can now estimate the energy M mat ter in the pure matter fields as a 
fraction of the total energy of the black hole in this two-fluid picture. We 
assume that up to the critical density p c the matter is pure and that the energy 
content contribution of the de Sitter fluid with an energy density p max starts 
to grow at p = p c . Then the matter content can be written as a two piece 
integral M maU er — 47r pr 2 dr + Air J^ c (3 (r) p c r 2 dr. We will assume (3 (r) to 
be simply linear in r with a form (3 (r) = ' which satisfies the requirements 
imposed on it above. Then substituting for (3 (r) and again noting that r 2 dr = 



1 2 
3 r r 9 



Mmatter — ^ 



we get 



+ 



l(r c + r )(r 2 c +r 2 )-h Q r -^ 



r 



(4.5) 



r 2 r Q In ( -B—Y 

u y V p max J 



From the definition of the density function we have that r = 

Since r — > r as p — > p max we have that for < r < r , r (p) is not well de- 
fined since r is now not single-valued in p. Further, we don't know the actual 
size of tq. As a result we shall assume, for the purpose of evaluating Eq. 4.5, 



that both r 



In 



and r — > ro as p — > p max hold so that 

fo (Pmax) — » 0. This is equivalent to having the de Sitter vacuum approached 
only as r — > 0. Thus the value of M mat ter calculated from here will be an up- 
per bound. Applying this on the integral J^ c (3 (r) p c r 2 dr part Eq. 4.5 implies 

\ (r c + r ) (r 2 + rg) - ±r (^f ) j^Pc so that Eq. 4.5 reduces to 



-^matter — 7^ 



^3 



3(a + l) 



(4.6) 
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Useofr c = 



-r%r g In 



and a = 2. 213 5 then gives = 0.479 17 [1 + 0.735 70] 
0.8316 - 0.83. 

This result agrees with the one found above in the quintessential picture. 



5 Conclusion 

In this paper we have suggested a possible model for a non-singular black hole 
as a product of gravitational collapse. It is an exact solution of the Einstein 
Equations, being Type [II, (II)] by Petrov classification. This model is based on 
a judicious choice we have made of the equation of state of the collapsed matter. 
At high densities this equation of state violates some of the energy conditions 
[2] [3] originally used to justify the existence of black hole singularities. In our 
model this violation leads to non-singular collapse. On the other hand, for 
the entire density parameter space the equation of state in our approach still 
satisfies the Weak Energy Condition p > 0, p + p r > 0, a basic requirement for 
physical fields. 

The solution depicts a spacetime with matter fluids in the outer layers (region 
II) which, as one moves deeper inside, give way to either a quintessential or a 
two-fluid region III. Region III, in turn, evolves into region IV, an inner-most 
core with de Sitter characteristics. The fluids in regions III and IV both provide 
the negative pressures needed to sustain the outer matter in static equilibrium. 
The solution has the required asymptotic forms, reducing to the Schwarzschild 
vacuum solution outside the matter fields and reducing to the de Sitter solution 
as one approaches the black hole center. The existence of region III renders 
the interface between matter fields and the de Sitter vacuum to join smoothly. 
This is because both the density p, the radial pressure p r and the tangential 
pressure p± profiles and their derivatives are continuous. This character of the 
fields makes the matching across associated interfaces natural. 

Using our model we have, in Section IV, offered two viable interpretations 
about the possible macro-state of the fields making up the total energy of 
the black hole. These interpretations are based on two pictures, namely a 
Quintessential Picture and a Two Fluid Picture. In both cases we have es- 
timated the fractional contribution of the regular matter-like fields as a fraction 
of the total black hole energy. We find for an upper bound the same value of 
~ 0.83 in both cases 3 . As a corollary, our model suggests that the lower bound 
for the amount of matter that is found in the "exotic" state can be as small 
as 17% of the entire configuration of the collapsed matter. Both the size and 
density of this de Sitter-like central core region will depend on the details of the 
microscopic physics that lead to the specific equation of state we have employed, 
which go beyond the scope of the present work. However, these considerations 
are consistent with our present notions of dynamical particle masses and sym- 
metry breaking, as the density of this state is much higher than those produced 

3 It is remarkable that the two different interpretations yield the same numerical values of 
the mass-fractions. 
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to-date in the laboratory. The precise value is uncertain but it has to be at least 
as high as m%, where tuh is the mass of the Higgs field, and it can be as high 
as the Planck density. Our entire non-singular configuration is well within the 
horizon of the black hole and in this respect the configuration is very different 
from some previous treatments, e.g. [6]. In these treatments, the field config- 
uration fills the entire volume interior to r — 2M with a fluid with p = —p, 
whose energy density for sufficiently large black holes can be smaller than that 
of water. This raises the difficult problem of converting matter from the usual 
equations of state to the "exotic" p = —p, under usual laboratory conditions. 

Finally, the model leaves some issues unresolved. First, the calculation for 
these mass fractions is based on assuming that in the region r < r < R both the 
pressure and the density are well-behaved, functions in the radial coordinate r, 
so that conversely r (p) = —r\r g In ) well defined in the entire parameter 

space of p. Since in our model p — > p max when r — > tq, then r (p) vanishes at tq 
forcing (in the mass calculation) the de Sitter state to appear virtually at the 
origin. It is in this sense that Mm ^" r — 0.83 is an upper-bound. One can reduce 

this fraction by using r (p) = —r\r g In {jr~ ) + r ° * n tnc calculations, instead. 
However since we don't know the radial extent the de Sitter field could fill at 
the core, or whether indeed it fills any space bigger than the Planck length we 
can not presently constrain 7V/1 "^" er from below. This is one issue for quantum 
gravity to settle. 

Secondly, our classical model can not distinguish between the two pictures 
offered to determine which choice of the matter macrostate is the correct one. 
However, the fact that in the region r$ < r < R the energy-momentum tensor 
elements p and p are well-behaved functions of the radial coordinate r may pro- 
vide an interesting insight. It suggests that one can associate a unique value of 
the energy-momentum tensor on each 2-surface in the family E = {£77 U £777}. 
One can further speculate that, such a configuration may even be quantizable. 
In such a case the information pertaining to the pre-collapse phase of the object 
would not be lost but would now reside on the (quantized) onion-like family of 
hyperfaces S inside the black hole. This is another issue for quantum gravity 
to settle. These results suggest a need for further investigation. 
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Figure 1: The quartic equation of state functional 
p (p) with m = 2, n = 1. At A when p — 0, then p = and = 0. The value of 
the parameter a — 2.213 (5) is chosen so that for the entire density parameter space 

p > the sound speed is not superluminal. The point B giving the maximum slope 

d 2 

corresponds to the maximum sound speed, = 0, chosen to correspond to the light 
speed, ^ = c 2 = 1. Point C is the critical point (p c ,p c ), and for p > p c , then p (p) 
is a decreasing function. At point D or p = p £ , the pressure temporarily vanishes.. 
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